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CONFLUENT HYPERGEOMETRIC ORTHOGONAL 
POLYNOMIALS RELATED TO THE RATIONAL QUANTUM 
CALOGERO SYSTEM WITH HARMONIC CONFINEMENT 



J. F. VAN DIEJEN 



Abstract. Two families (type A and type B) of confluent hypergeometric 
polynomials in several variables are studied. We describe the orthogonality 
properties, differential equations, and Pieri type recurrence formulas for these 
families. In the one-variable case, the polynomials in question reduce to the 
Hermite polynomials (type A) and the Laguerre polynomials (type B), respec- 
tively. The multivariable confluent hypergeometric families considered here 
may be used to diagonalize the rational quantum Calogero models with har- 
monic confinement (for the classical root systems) and are closely connected 
to the (symmetric) generalized spherical harmonics investigated by Dunkl. 



1. Introduction 

In this paper multivariable orthogonal polynomials are studied associated to the 
weight functions 
Type A (Hermite) 

(i.ia) a a (x)= n \ x 3-*k\ 2a ° n e ~^' 

1<j</e<71 1 <J <?l 

Type B (Laguerre) 

(1.1b) A B (x) = J] \(x j -x k )(x j +x k )\**> n \ Xj \^e-^l 

l<j<k<n l<j<n 

In the one- variable case (n = 1), the type A polynomials become Hermite poly- 
nomials (A A (x) = exp(— ujx 2 )) and the type B polynomials reduce to Laguerre 
polynomials of a quadratic argument (A B (x) = \x\ 2gi exp(— ujx 2 )). 

For both families we will exhibit (systems of) differential equations and Pieri type 
recurrence relations as well as the normalization constants that convert the polyno- 
mials into an orthonorma! system. Multiplication of the polynomials by the square 
root of the weight function, yields an orthogonal basis of eigenfunctions for the ra- 
tional quantum Calogero model with harmonic confinement and its generalizations 
associated to (the classical) root systems JCaj , pP[ . The connection between this 
orthogonal basis and the conventional (non-orthogonal) basis of eigenfunctions for 
the confined rational Calogero model (found by separating the quantum eigenvalue 
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problem in a 'radial' and a 'spherical' part) is described using the theory of Dunkl's 
generalized spherical harmonics with reflection group symmetry [ Dul , Du2| . 

The multivariable Hermite and Laguerre families associated to the weight func- 
tions A A 0) O^ ) and A B (x) ( jLlbj ) were introduced by Macdonald |M3| and 
Lassalle [ La2| , La3fl as a generalization (more accurately a deformation) of the pre- 
viously known special case in which the parameter go is being fixed at the value 
l/2gg,|Q. Recently, further insight regarding the properties of the poly- 



nomials considered by Macdonald and Lassalle was obtained in the context of a 
renewed study of the eigenvalue problem for the rational quantum Calogero model 
with harmonic confinement |D3|, UW, BF]. As it turns out, some of the results 



reported in the present work may also be obtained by combining results from pre- 
vious literature. For example, our evaluation formulas for the (squared) norms of 
the polynomials (cf. Theorem 2.2) can also be gleaned from [M3|, |La2|, |La3[ and 



BF], where expressions for these norms in a modified guise we re obtained by dif- 



ferent methods. (Specifically, if we make the norm formulas in [M3, |La2| , |La3|] and 



BF ] explicit with the aid of known evaluat ion formulas for the Jack symmetric 
functions at the identity due to Stanley |Sl], M4 |, then they are seen to be in cor- 
respondence with the expressions derived below in a completely different manner.) 
In all instances where overlap of this kind occurs (see the notes in Section 3.4), our 
approach provides an alternative, independent method of proof for the statements 
of interest. 

The paper is organized as follows. First, the confluent hypergeometric families 
associated to the weight functions A A (a;) (1.1a) and A B (x) ( LlbQ are defined in 
Section || and their main properties (orthogonality relations, orthonormalization 
constants, differential equations, and Pieri type recurrence relations) are formu- 
lated. In Section ^, we comment in some detail on the precise relation between our 
results and those obtained in previous literature. We will — in particular — take the 
opportunity to detail the connection between the multivariable Hermite/Laguerre 
families and the Calogero eigenfunctions as well as the relation to Dunkl's gener- 
alized spherical harmonics. In Section |] we provide the proofs for the statements 
in Section ^ by viewing the multivariable confluent hypergeometric Hermite and 
Laguerre families of the present work as a degeneration (viz. a limiting case) of 
certain families of multivariable hypergeometric orthogonal polynomials that were 
introduced in [D3] and investigated in more detail in [D4|. The multivariable hy- 
pergeometric polynomials relev ant to us here a re generalizations of the one- variable 
continuous Hahn polynomials [A2, AtSu , KS | (this corresponds to type A) and of 
the one- variable Wilson polynomials [|W, KS (this corresponds to type B) to the 
case of several variables. From a physical viewpoint, the multivariable hypergeo- 
metric polynomials in question are connected to the eigenfunctions of a difference 
(or 'relativisti c') c ounterpart of the rational Calogero models with harmonic con- 



finement |D2, D3, pi]. The transition from the hypergeometric to the confluent 
hypergeometric level corresponds to sending the difference step-size to zero. In 
this ('nonrelativistic') limit the difference Calogero model reduces to the ordinary 
Calogero model. Some technicalities needed to perform this transition at the level 
of the polynomials (which is established by controlling the convergence of the re- 
spective weight functions) are relegated to an appendix at the end of the paper. 



Note. Below we will always assume (unless explicitly stated otherwise) that the 
parameters go and <?i entering through the weight functions A A (x) (1.1a) and A B (x) 
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(1.1b) are nonnegative and, similarly, that the scale factor u> is positive. In principle 
it is possible to rescale the variables xi,... , s n so as to reduce to the case that u> is 
fixed at the value 1 (say). However, we have found it useful to keep the dependence 
on lo explicit in order to have a check on the scaling properties of our expressions 
and so as to suppress the emergence of numerical constants. 

2. MULTIVARIABLE HERMITE AND LAGUERRE POLYNOMIALS 

In this section the multivari able confluent hyp ergeo metric families associated 
to the weig ht functions A A (x) (|l.la|) and A B (a;) ( |l.lb|) are defined and the main 
properties of the polynomials are stated. The proof of these properties can be found 
in Section ||. 

2.1. Definition and orthogonality properties. Let m\, A G A denote the basis 
of symmetric monomials 



(2.1) 



m\(x) 



^6S„(A) 



A € A 



with 

(2.2) A = {A G Z" | Ax > A 2 > • • ■ > A„ > 0}. 

Here the summation in (|]l]) is over the orbit of A with respect to the action of the 
permutation group S n (which permutes the vector components Ai,... ,A„). We 
will also use the notation 

(2.3) m 2X (x)= ^ x\^---xf n , A e A 

M£S„(A) 

to indicate the basis of the sym metr ic mo nom ials that are even in the variables 
xi, . . . ,x n . The monomial bases (2.1) and (2.3) inherit a partial ordering from the 
dominance type partial ordering of the cone A (2.2) that is defined for A,/i 6 A by 

A < (j, iff > A, < > n-j for k = 1, 



(2.4) 



J2 x j< ^ 

l<j<k l<j<k 



(A < jJ, iff A < /x and A^/i). 

Let ( •■ -)a and (•, -)b de note the L 2 inner products over 



A A (x) ( 1.1a ) and A B (x) ( 1.1b ), respectively. So, explicitly we have 
(2.5a) 



with weight function 



(for f,g in L 2 ( 
(2.5b) 



(/.P)a = 
1", A A dxi 

(f,g)B = 



f{x) g(x) A A (x) dxt 



■ dx n ) ) and 



f(x) g(x) A (x) dx\ ■ ■ ■ dx n 



(for f,g in L 2 (M. n , A B dxi ■ ■ ■ dx n )). After these notational preliminaries, we are 
now in the position to define the multivariable confluent hypergeometric families 
associated to the weight functions A A (x) ( 1.1a ) and A B (a;) ( |l.lb| ). 

Definition. The type A (or multivariable Hermite) polynomials p A (x), A G A are 
the polynomials determined (uniquely) by the conditions 

A.l pf(x) = m\(x) + c a,m to m( x ): 



/iGA,/j<A 



A.2 (jp A ,m fJ ,) A = if ^<A. 



4 



J. F. VAN DIEJEN 



Similarly, the type B (or multivariable Laguerre) polynomials pf (x) , A G A are the 
polynomials determined (uniquely) by the conditions 

B.l pl{x)=m 2X {x)+ c i» m z»( x )> c f,M eC 

B.2 (pf,m 2fl )B = if n < A. 

The type A polynomials p^ (x) , A e A constitute a basis for the space of permu- 
tation invariant polynomials in the variables xi, . . . ,x n and the type B polynomials 
pf(x), A G A form a basis for the even subsector of this space (i.e., the subspace of 
symmetric polynomials in x\ , . . . , x%). The following theorem says that the bases in 
quest ion are orthogonal with respect to the inner products (•, -)a ( 2.5a ) and (•, -)b 
( [2.5b|) , respectively. 

Theorem 2.1 (Orthogonality). Let A,/x G A (2.2). We have 



(Px > P°)c= ■■ p%(x)pC(x)A c {x) d Xl ---dx n 

J OO J — CO 

= if A^/x, 

where C stands for A or B. 

For weight vectors A and /i that are comparable with respect to the partial order 
( p^ ) the orthogonality of p^ and p^ follows immediately from the definition of the 
polynomials. Theorem 2.1 states that the orthogonality relations in fact hold for 



general weight ve ctor s A, fj, G A (|2.2| ) (not necessarily comparable with respect to 
the partial order ( |2.4|) ). 

In order to orthonormalize the bases {P\}\eA arid {p^}a<eA, it is needed to 
evaluate the integrals for the (squared) norms of the polynomials. 



Theorem 2.2 (Norm formulas). Let A G A (|2.2|). We have 

/oo />oo 
•••/ \px{x)\ 2 A A (x) dxi ■ ■ ■ dx n 
-oo J — OO 

_ (27r)"/ 2 n! 

(2w)l A l+fo«(»-l)/2+n/2 

T((fc -j + l) 5o + Aj - A fe ) T(l + {k-j- l) 5o + Aj - A fe ) 



x 



n 



, . , A „ T((fc - j)g + X 3 - Afe) T(l + (k - j)go + A, - \ k ) 
x H r(l + (n-j)go + Xj) 

l<j<n 

and 

/oo />oo 
•••/ \p*(x)\ 2 A B (x)d Xl ---dx n 
-oo J —oo 

n! 



w 2|A|+gon(n-l)+(ffi+l/2)n 

r((fc - j + l) 5o + A, - Afe) T(l + {k-j- l) 5o + A, - A fe ) 



x 



n 



, ( , „ r((k-j)g + X 3 -Afe)r(l + (/c-j).g + A J -Afe) 



x J] r(l + (n-j).9o + A J )r((n-j)<?o + .gi + l/2 + A J ) 
(where |A| = Ai + •• • + A„ and T(-) denotes the gamma function). 
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For A = the polynomials p A and £? reduce to the unit polynomial (p A {x) 
Po(x) = 1). The formulas in Theorem 2.2 are then seen to simplify to 



(2.7a) 
and 

(2.7b) (1,1) B 



/oo poo 
■■■ A A {x)dx!---dx 
-oo J — oo 

(2tt)"/ 2 



TT r(l+Jff ) 
(2 W )90«(n-l)/2+n/2 11 T(l + ) 



A B (x) dici ■ • • dx r , 



1 



w 9o"(n-l)+(ffi+l/2)n 



n 



l<j"<n 



r((j-i)go + gi + i/2)r(i + j go ) 
r(i + 5o) 



The integra ls in ( 2.7a ) and ( 2.7b ) amount to integrals evaluated by Mehta and 
Macdonald [Me, Ml|. More precisely, Mehta conjectured the closed expression for 
the value of the integral (1, 1)a and Macdonald generalized this conjecture (in terms 
of root systems) therewith also including the integral (1, 1)b- Both the integration 



formulas for (1,1) a an d (1, 1)b m ( 2.7a ) and ( 2.7b ) were then proven in [Ml] by 
viewing them as limiting cases of an integration formula due to Selberg 
also the introduction of | Al[). 



M ( cf - 



2.2. Differential equations. First we introduce two families of commuting dif- 
ference operators D A g , . . . ,D A g and Df a, ■ ■ ■ >^n,8- The tv P e ^ operators are 
given by 



(2.8a) 



D 



E 



U J a + r\J±,r-\J + \-\J-\ v J+J-;J^nJt e 



J+>J_C{1,... ,r>} 



J +n j_=0, |J+| + |J_|<r 

r = 1. . . . , n, with 



-(d J+ -dj_) 



V A 

V J + ,J-;K 



II .<**' II 

jeJ+ jeJ- 

n wA (^) n wA (- x i) n yA ^ - ^> A o^ - ^ - 

jeJ+ je.J- je.J+,j'£J- 

x Y[ v A (xj - x k ) Y[ v A (x k - xj), 



keK 



keK 



n A 

K,p 



lEL + leL- 

Yl v A {xi - x v )v A {x v ~xi+ i/3) 



|L + | + |L_|=p 



VeL- 



Yl v A {xi-x k ) Yl v A {x k -xi) 



keK\L + UL. 



ieL_ 
keK\L + UL. 
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and 

flgo 



v A ( Z ) = fi+^y w A { Z ) = 

The type B operators read 

(2-8b) D B ,= £ C^^i^e^ 



JC{1,... ,n},0<|J|<r 
£j=±l, JGJ 



r = 1, . . . , n, with 

= IJe 6 '^ foe {-1,1}), 



e ' 



Kj,k = Yl^(sjXj) Y[ v B (s j x j +£ j >x j >)v B (e j x j + £ f x f -i0) 

j<j' 

x JJ v B (ejXj + x k )v B (ejXj - Xk), 



keK 



U B 



(-l) p Yl (j[w B (eixi) ]]_ v B (eix l +e v xi,)v B (-eix l -ei l xi, +if3) 



LCK,\L\=p IGL l,l'€L 

£I=±1, l£L 1<V 

x Yl v B (em + x k )v B {eixi - x k ) 

LeL 
k£K\L 



and 



^ z)= ( 1 + ^l\ w B (z)=(l + (1 + i/3uz) 

\ IZ ) \ IZ / 



In ( 2.8a ) the sum is meant over all disjoint pair s of i ndex sets J+, J_ G {!,■■■ , fi} 



with the sum of the cardinalities being < r. In ( 2.8b ) the sum is over all index sets 
J C {1, ... , n } of cardinality < r and all configurations of signs £j € { + 1, — 1}, 
j € J. We have furthermore assumed the conventions that empty products are 
equal to one and that and U B are equal to one for p = 0. The exponentials 

exp(±jg|-) act on analytic functions of X\,.. . ,x n by means of a shift of the jih 
argument in the (for (3 real) imaginary direction 

±(- — ) 

(e ; ax J /)(xi, . . . ,x n ) = f(xi, . . . ,Xj-i,Xj =F i(3,Xj+i, ■ ■ ■ ,x n ). 



Thus, the operators D A g ( 2.8a ) and D B ( 2.8b ) are analytic difference operators 



of order 2r in exp(|^-), . . . ,exp(f ^f-). 

If we act with D^^ (C = A, B) on an (arbitrary) analytic function / of the 
variables X\, . . . ,x n , then we end up with an expression that is holomorphic in the 
step size parameter f3 around (3 = 0. The coefficients of the Taylor expansion in 
(3 around zero can be written in terms of partial differential operators applied to 
the function /. We will call the first nonzero differential operator in this expansion 
the leading differential part of the difference operator. One obtains the leading 
differential part D^ of D®n (C — A, B) by expanding the analytic difference oper- 
ator in powers of (3 using the formal identity exp(±jg|-) = $Z^_o(±^ a§~) m / TO " 
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The result is a formal expansion of the difference operator in terms of differential 
operators that gets its precise meaning (as a Taylor expansion) upon acting with 
both sides on an (arbitrary) analytic function. 

The following theorem describes the structure of the highest order symbol of 
the leading differential parts D A Q , . . . , D A and Df , . . . , D^ for the difference 



operators D A ( 2.8a ) and ( 2.8b ), and states that the differential operators of 
interest are diagonal with respect to the bases {P\}\eK and {p\}\eA, respectively. 

Theorem 2.3 (Differential equations), i. The formal expansion of the difference 
operator D^a (r = 1, . . . ,n; C = A, B) in the step size parameter (3 has the form 

D% = B%^ +0(f3 2r+1 ) 

with 

^o = (-ir E n£ 

\J\=r 

(where l.o. stands for the terms of lower order in the derivatives d/dxi, . . . , d/dx n ). 

ii. The leading differential parts D± Q , . . . , D~j commute and are simultaneously 
diagonalized by the multivariate Hermite ( type A ) and Laguerre ( type B ) polyno- 
mials 

BC pC = E r c (\)pC, AeA, 
with the corresponding eigenvalues given explicitly by 

E^{X) = (2 U y J2 H X > and E*(X) = (^y J2 U X r 

JC{1,... ,n} j£j JC{1,... ,n} j£j 

\J\=r \J\=r 

Corollary 2.4. We have that 

lim (3- 2r (x) = E r c (X) p°{x) , X e A. 

Corollary 2.5 (Symmetry). The leading differential operators -Df i • • • i ^nO ma P 
the space of permutation-invariant polynomials in x\, . . . ,x n (type A) or x\, 



(type B ) into itself and are symmetric with respect to the inner product (•, -)c (2.5a), 

( B) 

(£ r c mf , mfic = (mf , D% m C) c , X, /* € A 

(with m^(x) = m\(x) (2.1) and = m2\(x) (2.3) ). 



The eigenvalue equations (L>^ p£)(x) = lim/3 2r (D^ p Px)(x) = E?(\)p%(x), 

r = 1, . . . , n, constitute a system of differential equations for the polynomials (x) , 
AeA. It is instructive to exhibit these differential equations in a more explicit 
form for the case r = 1 (which amounts to the order of the equation being equal to 
two). For r = 1, the difference operators D A p ( |2.8a ) and Dfp ( [2.8b]) reduce to 



(2.9a) D A = (™ A (xj) II ^( Xj - x k )(e* I - 1) 



l<j<n l<k<n,k^j 

w A (-x 3 ) 11 v A (-x j + x k )(e~ T ^ - 1) 
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and 



(2.9b) D B 



E 

l<j'<n 



v B (xj - x k )v B (xj + Xfc)(e 1 *i - 1) 



w a (~ Xj ) [[ v a {-x j +x k )v a {-x j -x k ){e Xs i -1) 

l<fe<n,fc^j 



(with v A (-), w A (-) and « B (-); U,B (') taken the same as in ( |2.8a| ) and ( 2.8b| ), re- 
spectively). A formal expansion in (3 reveals that „ = D^ 1 + 0(/3 3 ) with 



(2.10a) D A = (-7^ + 2uJX ^)- 2 ^ E 



l<j<n 



(2.10b) D B fi = J2 

l<j<n 



dx 3 

dxj 



25] 



dx 3 
1 d 



2ujx 



l<j<k<n 

d 



-(- 

Xk V dx-i 



d 
dx k 



-2.9o E 

l<j<k<n 



.( 



3 u - h 3 
8 



dx q 



x j x V 



d 

dxk 



) +— ( 

/ i, + Xk V 



Xj + Xk \ dxj dxk 



Formula (2.10a) and (2.10b) combined with the corresponding expressions for the 



eigenvalues taken from Theorem 2.3 

(2.11) £ A (A) = 2w(A x + • • • + A„), E B {\) = 4w(A 1 + --- + A„), 

render the second order differential equations Z? A p A = _E A (A)p A and Df p B = 
Ef(X)p B in a fully explicit form. 

In principle it is possible to determine for given r the differential operator D^ 
algorithmically as the leading part of the (explicit) difference operator -D^s ( 2-8a ) 
,(2. 8b) (via the formal expansion in the step size parameter 0). It seems a rather 
nontrivial combinatorial exercise, however, to derive along these lines a closed ex- 
pression for the differential operator D^ Q for general r. 

2.3. Pieri type recurrence formulas. To describe the recurrence relations for 
the multivariable Hermite and Laguerre polynomials it is convenient to pass from 
monic polynomials to a different normalization. Let 

P A (z) = cM(z), P B (x) = c B p B (x) 

[{k- j)go]x J -x k 



(2.12) 




with 




(2.13a) 




(2.13b) 


of 



n 



l<j<k<n 

(- W )W 



[(1 + 

n 



fc - 3)go}k :i -x k ' 

p-i)5o]A,-A k 



[(1 



l<j<k<n 

x TT 

ii<„ [(n_j) - 9o+5i + 1/2] ^ 



i 



where we have employed the Pochhammer symbol defined by [a]o = 1 and [a]i 
a(a + 1) • • • (a + I - lj for I = 1, 2, 3, ... . 
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Theorem 2.6 (Normalization). The normalization of the polynomial P^ (x) ( [2.12 ) 
is such that 

lim a~ |A| P^(al) = 1 [with 1 = (1, . . . , 1)) 

a — >oc 

and 

P A B (0) = 1 {with 0= (0,... ,0)). 



The next theorem describes an expansion formula for the product of Pjp (x) ( [2.12 ) 
and the first elementary symmetric function in x\, . . . ,x n (type A) or in x\ , 



(type B). Formulas of this type are often referred to as Pieri formulas [M4, |St[. 
(More generally, Pieri formulas are relations in a commutative algebra describing 
the expansion (in terms of a basis) of products between the basis elements and a 
set of generators for the algebra.) For n = 1 the formulas in the theorem reduce to 
classical three-term recurrence relations for the one- variable Hermite and Laguerre 



polynomials (cf. Section 3.1) 



Theorem 2.7 (Pieri formulas: simplest case). 

The (renormalized) multivariable Hermite and Laguerre polynomials P^(x) ( |2.12 ) 



satisfy the recurrence relations (ej denotes the jth unit vector in the standard basis 
ofWL n ) 



E 4tf(*) = 

l<j<n 

E {vfP> B +ej (x) - (Vf + V»)P x B (x) + \ , («)) 



with 



l<j<n 

a TT fi i 9o 



% A = n ( i+ 

l<k<n,k=£j 



{k - j)go + Xj - > 
{n - j)g + Aj i r / go 



-3 0„, 11 



l<fe<n,K7=j J 



Vf = ((n — j)go + gi + 1/2 + Aj) JJ (l 

i<fc<?i,fc^j 



,9o 



( fc - j).9o + Aj - A fe 



,9o 



l<fc<n,fc^j v /"u 

One word of caution is at place here. It may of course happen that for certain 
A £ A and j € {!,••• , the vector A + ej (or A — ej) does not lie in the cone A 



(2.2). For such boundary situations the polynomial PS_ e (or Pjp_ e .) is not defined 
and it might a priori seem that r.h.s. of the recurrence relation does not make 
sense in this case. It is not difficult to verify, however, that in these situations the 
coefficient Vp (or VSj) in front of Pf +e . (or Pf_ e .) vanishes. (Indeed, for A € A 
we have that A + e j A if Aj_i = Aj and that A — ej ^ A if Aj = Aj+i or if 
j = n and X n = 0. In the former case we pick up a zero in from the factor 
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1 + 9o((k — 3) 9a + -\j ; — Xk) 1 with k = j — 1 and in the latter case we have a zero 
in V_j from the factor 1 — 9o((k — j)go + Xj — Xk)^ 1 with k — j + 1 or from the 
factor (n — j)go + Xj with j = n, respectively.) 



The next step is to generalize the expansion formulas of Theorem 2.7 to (multi- 
plication by) arbitrary elementary symmetric functions. To this end we introduce 
(for r = 1, . . . ,n) 

(2.14a) E^(x) = J2 II^> 

JC{1 ™}j£J 

\J\=r 

(2.14b) Ef{x) = {- u ) r II 



T 2 



JC{1,... ,n}jeJ 
\J\=r 

It is clear that the products £y (x)Pjp(x) can be written as a linear combination 
of Pp(x) with /.t < A + e\ + ■ ■ ■ + e r (this is immediate from the structure of 
the monomial expansion of P\{x) and the fact that such expansion formulas for 
these products evidently hold if we replace the polynomials Pjp (x) by their leading 
monomials m\(x) (C = A) and n%2\(x) [C = B)). It turns out that many of the 
coefficients in the expansion E^P^ — X^<a c m-^/? are m ^ &c ^ zer0 (f° r r = 1 this 



is of course apparent from Theorem 2.7). The following theorem provides detailed 
information on the structure of the terms entering the Pieri type expansion of the 
product between the basis element P\(x) and an arbitrary elementary symmetric 
function E°(x) fl2.14a| ), (|2.14b| ). 



Theorem 2.8 (Pieri formulas: general structure and leading coefficients). 

The (renormalized) multivariable Hermite and Laguerre polynomials P^(x) ( p. 12 ) 

satisfy a system of recurrence relations of the form (ej = J2jeJ e i) 

E r c (x)pC(x) = 

E Wf +tJ _. r (X) P? +ej+ . ej _ (x) r = 1, . . . ,n. 

J + ,J_C{1,... ,»}, A + e J+ -ej_ £A 

J + nJ_=0, |J + | + |./_|<r 

TTie coefficients Wj + j . r (X) that correspond to index sets J + , J_ with the sum of 
the cardinalities \J-\-\ + \J—\ being equal to r are explicitly given by 

Wj C +>J _. <r (X) = vf +tJ _. {J+UJ _ )c 

where 

yA -Q f( n ~ 3)90 + Xj • 

jeJ- 
x 



5o \ tt / 1 5o 



+ (i' - i)5o + Aj - A, 



TT 1 + 29 TT 1 

1 1 V (k-j)g + Xj - X k J j£_ V (k - j)go + X 3 - X k 



je.J 

k£K keK 
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and 



V 



J + ,J-;K 



II (( n - 3)90 +9i + 1/2 + Xj) I] ((" - 3)30 

9o 



A,) 



>< n 

jeJ+,j'eJ 



x 



n 

keK 



1 



1 + Tl— •> 

9o 

(k - j)go + Xj - x k 



Xf 



n 

keK 



1 



.9o 



1 + (f - ^').9o + Aj 

5o 

(k - j)go + X 3 - X k 



(with the convention that empty products are equal to one). 



Theorem 2.8 constitutes a partial generalization of Theorem 2.7. For r = 1 
the structure described in Theorem |2.8| is compatible with that of the formulas 
in Theorem 2.7 and we furthermore recover the coefficient in the r.h.s. of Pjp +e . 



(viz. Vp 



'{j}Mh...,n}\{j} J and P ^ ( V1Z - r ®,{jh{i,-ri\{o} J uuu 11UU ^ U1 ^ 
(which — according to Theorem |2.7| — happens to be zero for C — A and equal to 



V 0,{j};{l,...,n}\{j} 



) but not that of Pjp 



jV{i},0;{l,...,n}\{j} 



*}\{j} 



) for C = B). 



Even though Theorem |2.8| is not completely explicit, as it does not tell us the 
expansion coefficients Wj + j_. r (X) for | J + | + |J_| < r, it is still useful in its present 
form. For instance, the theorem implies (together with the orthogonality) that 

(2.15) (4 c P A c ;P c )c = 



V J+,J_;(J + UJ_) C 



if /i 7^ A + ej + 
if /i = A + ej + 



with | J + 1 
with |J_|_| 



|J_|<r 
\J-\ = r 



(where J + , J_ 6 {1, 
identity (£?°Pf ,P!" 



A+e, 



, n} such that J+ n J- = 0). When applying (2.15) to the 

>c one arrives at a system of 



C nC 



A+e, 



recurrence relations for the squared norm of Pp 



(2.16) 



.,r},0;{r,...,n}(A) (P\+e {1 , 



,P 



A+e, 



{i, 



,r};{n...,n}( A + e fl-.,r}) (P\ , P\) 



(r = 1, . . . , n). The recurrence relations in ( 2.1 6| ) determine (Pjp, Pjp)c uniquely in 
terms of (1, l)c (because the (fundamental weight) vectors lT -}> r — 1, ■ ■ ■ , n 
positively generate the cone A (2.2) and the coefficient r y g.r r n y(X) 7^ for 

A € A). This observation gives rise to an alternative (constructive) proof of the 



Indeed, by using the property c 
( fU6| ) in the monic form 



Vfe { i t . 



norm formulas in Theorem 2.2 different from the proof presented in Section 4.1 

ji(A) one rewrites 



"},0i{r,. 



(2.17) (pj+e {li ... ir} ,^ + 



e {l.- 



. ,r}' 

,r},0;{r,... ,7i}(A)V0,{l, ..., r };{r,. 

which upon iteration and matching of th e i nitial conditions s o as to reduce for A = 
to the Mehta-Macdonald formulas ( 2.7a ), ( 2. 71 ) (cf. Section 4.1) leads to the norm 
formulas of Theorem 2.2. 

The last theorem (below) provides a complete (explicit) description of the ex- 
pansion coefficients Wj j_. r (A) for the type B case (thus including the coefficients 



l }(A + e {li ... )r} )(p A ,p A >c 
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corresponding to index sets with \J+ \ + |J_| < r). This renders the system of re- 
currence relations of the type given by Theorem 2.8 in a fully explicit form for the 
multivariable Laguerre family. 

Theorem 2.9 (Pieri formulas: explicit expansion coefficients Laguerre cas e) . 
The coefficients in the recurrence relations of the type described b y Th eorem 2.8 are 
for the renormalized multivariable Laguerre polynomials P^{x) (2.12) given by 



with V 



W J + ,J-;rW — V f+,J-;(J + UJ-)" ^(J + UJ-)",r-\J+\-\J-\ 

a taken the same as in Theorem \2. q and 



J+,J-\ (J+UJ_) 

^!, P = (-i) p x 



E ( II « n - 0.9o + .91 + 1/2 + Xj) J] ((n - l)g + X/) 

k . Li nL _ =0 \IgL-l- /GL_ 



L + ,L__ CK, 

\L+\ + \L_\=p 



X 



n 



.90 



x n(i+ 



(l> - 0.9o + A, - A, 

A n(i 



i 



,9o 



1 + {V - l)g + X t - V 
3o X 



IEL + 



(k - I) g + X[ - X k 
( with the convention that = 1 for p = 0) 



(k - 0.9o + X[ - X 



r) 



3. Comments 

3.1. The special case n = 1. In the case of one single variable (n — 1), the weight 
functions reduce to 

(3.1) A A (x) = e-" x \ A B {x)^\x\ 29l e- LOx2 . 

The polynomials then become monic Hermite polynomials (type A) and monic 
Laguerre polynomials of a quadratic argument (type B), which can be written 
explicitly in terms of a terminating confluent hypergeometric series [A.bSt| 



(3.2a) p$(x) 
(3.2b) p B x (x) 




for A even 
for A odd 



(with 1F1 ( " ;z) = V [u,Jm and A = 0,1,2,...). The norm formulas, 



Jo|. 

differential equations and recurrence relations reduce in this special situation to 
classical formulas for the one- variable Hermite and the Laguerre polynomials (notice 
however that the scale parameter u> is usually taken to be equal to one and that 
our normalization differs from the standard one) 
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Norm formulas (cf. Theorem 



(3.3a) 
(3.3b) 



\pUx)\ 2 A B (x)dx = 



/OO 
\p*(x)\ 2 A*(x)dx = 
-OO 



XU 



2 A w A+l/2 



A!r( 5l + l/2 + A) 



u) 2X+g 1 + l/2 

Differential equations (cf. Theorem 2.5 and Eqs. ( p. 10a ). (2.10b) and (2.11)) 



(3.4a) 

(3 ' 4b) W x dx 

Recurrence relations (cf. Theorem |2.7| 

A 



d 2 B 2gx d B 



d o 
dx FX 



4wAp B 



(3.5a) 



xPt = 



P 



A+l 



2uj 



-P 



A-l 



(3.5b) -cux 2 Pj? 



(ffi + 1/2 + A)P A B +1 - ( 9l + 1/2 + 2A)P A B + APf_! 



with P£(x) 



Px(x) and P B (a:) = rJ^r£mn; P\ ( x )' ^ ne norma lization properties 
for P A (x) and Pjf (x) in Theorem |2.6| state that limo.^00 a _A P A (a) — 1 and that 
P? (0) = 1 . In the present situation these properties are immediate from the explicit 
confluent hypergeometric representations in (3.2a) and (3.2b). 

3.2. Eigenfunctions for the rational Calogero model in an harmonic well. 

Most of the results stated in Section admit an interpretation in terms of certain 
exactly solvable quantum mechanical n-particle models on the line. Specifically, 
conjugation with the square root of the weight function A c (x) (1.1a), fll.lbD trans- 
forms the second order differential operators Dp ( p. 10a ), (2.10b) into Hamiltonians 
for the rational quantum Calogero models with harmonic confinement associated 
to the classical root systems Ca, OP] 

(3.6a) ff A = (A A )^L»f (A A )-^ = 



y (-— 

4r*. V dx 2 



2 2 



(3.6b) 



l<j<n 

Hf = (A B )3^f(A B ) 

d 2 



2^0(30-1) y (xj-xkY 

l<j<k<n 



-Ei 



y 

l<j<n 



-2.9o (.90 - 1 



y 

l<j<k<n 

with Eq = wn(l + go(n — 1)) and Eq 
Theorem 2J3 that the functions 

^(x) = (A c (x))^ p°(x), AG A 



(3.7) 



( a _ 2 , ,\ 

XTdx 2 . +91(91-^, +uj Xj j 

((xj - x k y 2 + (xj - x k )~ 2 

= wn(l + 2go(n — 1) + 2gi). It is clear from 
(C = A,B) 



constitute a basis of eigenfunctions for (3.6a), (3.6b) with the corresponding 
eigenvalues given by E^{\) ( p. 11 ). In its present form these eigenfunctions for the 
rational Calogero model with harmonic term were introduced in |D3| , BF] and also 
(for type A) in [UW]. The orthogonality (in i 2 (R", dx\, . . . ,dx n )) of the basis 
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ipx( x ) (3-7) follows from Theorem 2.1 and the orthonormalization constants can be 



read-off from Theorem 2.2 



Historically, the study of the eigenvalue problem for the type A Hamiltonian 
(3.6a) was initiated by Calogero, who computed the spectrum and determined 
the structure of the corresponding eigenfunctions to be a product of the ground- 
state wave function ip A (x) — (A A (x)) 1 ^ 2 and certain symmetric polynomials in 
Xi,...,x n [Ca|. To be precise, Calogero considered a translationally symmetric 
n-particle system with a potential of the form 



V(x) 



E 

l<j<k<n 



(G (Xj 



Xk) 



Xk) 



which is seen to be equivalent to the type A system above up to a simple center of 
mass motion {J2j <k {xj — x k ) 2 — n J2j x 2 — (J2j x j) 2 )- Explicit expressions for the 
Calogero eigenfunctions were found for small particle number (n < 5) by Perelomov 
and Gambardella []P(J |G|. More recently, it was observed |BHV, BHKV| that for 
arbitrary particle number n it is possible to construct a basis of eigenfunctions for 



H A (3.6a) with the aid of certain creation and annihilation operators and a~ 
(j = 1, . . . , n) t hat are built from the differential-reflection operators introduced 
by Dunkl |Du3|. In a nutshell: the relevant operators are obtained starting 



from the usual creation/annihilation operators (Tgf - + &Xj) for a system of un- 
coupled (bosonic) harmonic oscillators (this corresponds to go = 0) by replacing 
the partial derivatives by the corresponding Dunkl differential-reflection operators 
(associated to the root system A n ^i). If one acts on the ground-state wave function 
^o(x) = (A A (x)) 1 ^ 2 (which is annihilated by a^, . . . , a~) with an arbitrary sym- 
metric homogeneous polynomial of degree I in the creation operators af , . . . , , 



then one winds up with a symmetric eigenfunction of H A (3.6a) with eigenvalue 
21lu. By taking for the symmetric polynomial in question the Jack symmetric func- 
tion J\(x; 1/ 'go) ||St , M4], one arrives (up to a normalization factor) precisely at the 
eigenfunction V A (x) (O fUW]. 



The basis of eigenfunctions of the form ^^(x) (3.7) for H± (3.6a), (3.6b) is also 
very special in that it simultaneously diagonalizes the higher-order quantum inte- 
grals for the Calogero model. Specifically, by applying the similarity transformation 
Di — * Hi also to the higher-order differential operators in Theorem 2.3, one 



obtains a complete set of commuting quantum integrals for the Calogero model of 
the form 



(3.8) 



H? = (A C )?D A (A C ) 
(-1) ^ 



Je{i, 



\J\-- 



n— 



+ l.o. 



(r = 1, . . . ,n). 



It is immediate from Theorem |2.3| that the functions ^^(x) (3.7) constitute a basis 
of joint eigenfunctions for the differential operators H^ (|3~g ), r = 1, . . . ,n. Since 
the corresponding eigenvalues -Ef(A),... ,E^(X) separate the points of the cone 
A ( |2.2| ), it follows that this property actually determines such basis uniquely. In 
other words, ambiguities in the possible choices for the basis of eigenfunctions of fff 
(3.6a), ( |3. 6b ) caused by the degeneracy of the spectrum (2.11) get eliminated by re- 
quiring the eigenfunctions to be a basis of joint eigenfunctions for . . . , H^ . The 
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other important (and restrictive) property of the basis ^>^(x) (3.7), viz. its orthog- 



onality, implies — together with the real-valuedness of the eigenvalues E® (A) — that 
the differential operators are (essentially) self-adjoint (since unitarily equivalent 
to real multiplication operators) in the Hilbert space of permutation-invariant (type 
A) or permutation- invariant and even (type B) functions in L 2 (R™, dx\ ■ ■ ■ dx n ) (cf. 
Corollary U). 

One would like to cast the differential operators for arbitrary re {1, . . . , n} 
in a more explicit form. Explicit formulas for n-independent commuting differential 
operators generating the same algebra as our operators H B , . . . , H B appear as spe- 
cial cases of the commuting families presented in [OOS, OS . In order to determine 
the relation between the operators H B ( |3.8| ) and the type B differential operators 
in | OOS, OS] explicitly (for all r), one would have to know the eigenvalues of the 
latter operators on the basis ^ B {x) ( |3.7| ). For the type A we are not aware that 
similar explicit formulas for a set of commuting differential operators generating 
the same (commuting) algebra as H A , . . . , H A have been reported in the literature 
(except for oj — | OOS| ). However, it is known that such commuting differential 
operators may be characterized in terms of Dunkl's differential-reflection opera- 



tors [Po| and recently the eigenvalues of the thus obtained differential operators on 

(thus allowing a comparison with our 



the basis $£(x) @ 

were determined 
differential operators H A ) 



3.3. Relation to Dunkl's generalized spherical harmonics. An alternative 
approach towards the solution of the eigenvalue problem for the second order dif- 
ferential operators D± ( 2.10a ), ( 2.101: ) is to separate the eigenvalue equation in a 
'radial' and a 'spherical' component. In essence this is the method used by Calogero 
|Ca[ to obtain the eigenfunctions for the type A Hamiltonian H A (3.6a) (cf. the 
comments in the previous subsection). Specifically, if one substitutes an Ansatz 
function of the form R c (r)Yi(x), where R c (r) is a function of r = \/ x\ + ■ ■ ■ + x\ 
and Y l (x) is a permutation symmetric homogeneous polynomial of degree I in 

X\,.., 



(C = A) or x\,. 



(C 



B), then it is seen that this yields an eigen- 
function of (|2.10aD , fl2.10b| ) with eigenvalue E c if R c (r) and Yp(x) satisfy 



(3.9a) 
(3.9b) 



d 2 R A 

dr 2 
d 2 R B 
dr 2 
and 

(3.10) 

with 

(3.11a) L A = 
(3.11b) L E 



,^ = (E A -2ul)R A , 
r J dr 

41 + g n(n - 1) + 2ng 1 + n - 1\ dR B 



^ r _ 2l+g Q n(n- 1) + n - 1 ^ 
(2lut 



dr 



l c yF = o 



V- 9 2 „ x , 

E d^ + 2 9o E 



(C = A,B) 



1 



d 2 



dx 2+2gi 




/ d d \ 
Xj — Xk V dxj dxk ' ' 



= {E B - 4^2)^ 



E v 

l<j<n \ 3 3 

1/9 d \ 1/9 
%j — Xk \dxj dxk' Xj + Xk \dxj dxk- 



+250 J2 

l<j<k<n 



d 
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The 'radial' equations (3.9a) and (3.9t ) are confluent hypergeometric type equations 
that admit polynomial solutions for E A = 2w(l + 2m) and E B = Alo(1 + m) (with 
in G N) given by Laguerre polynomials in r 2 (cf. Eqs. (3.2b) and ( 3.4b[ )) 



(3.12a) Ri(r) 



(3.12b) R B m {r) 



[l + (n-l)(l+n 5 o)/2]r 



x iiq 



1 1 I + (n-l)(l + nso)/2 
[21 + n(l/2 + gi + (n - l)go)] m 



— m 

2l + n(l/2 + gi + (n- l)g ) 



(Here we have chosen the no rmal ization chosen such that i?^ (r) is monic.) 

The 'spherical' equation ( 3.10| ) was studied (for type A) by Calogero |Ca| and 
in furt her d etail and more generality (therewith also including the type B) by 
Dunkl |Dul|. Let Vf be the space of homogeneous symmetric polynomials of 



degree ! in x\, . . . , x n (C — A) or xf, . 



(C = B) and let U\ C 



Vp be the 

subspace of polynomials satisfying ( 3.1 0| ) . The polynomials in Tlf are referred to 
as (symmetric) generalized spherical harmonics. For go,gi = these generalized 
spherical harmonics reduce to ordinary harmonic polynomials in M™. It follows from 
Dunkl's theory in |Dul (see also [ Du2{ for the extension to the nonsymmetric case) 
that dim(Hf) = dim(7>, A ) - dim("P^7) and that dimfftf ) = dim("P B ) - dim("P ; B 1 ). 
The upshot is that each polynomial p^ may be written uniquely in the form 



PI/2] 



(3.13a) p£(x) = R ™( r ) Y \\\-2m(x), 



Y 



|A|-2r. 



e H 



|A|-2r, 



m— 



(3.13b) p B x (x) = J2 R *W Y W-m(x), 



r |A|-m fe n \\\-m 



m=0 



(where [•] represents the function that extracts the integer part). Indeed, the func- 
tions of the form in the r.h.s. of ( 3.13a| ) and (3.13b) are eigenfunctions of 
( [2.10a| ), (2.10b) corresponding to the eigenvalue 2ui\X\ (type A) and 4w|A| (type 
B), respectively. Furthermore, the functions of this form span a space of dimension 



y 



|A|/2] 
=0 



dim(Hf A| _ 2m ) = dim(pA) and £ 



, B ,), which is 



precisely the multiplicity of the eigenvalues -Bp (A) (2.11). 

The formulas ( |3.13a ), ( 3.13b| ) describe the relation between the Calogero type 
eigenfunctions of the form R c (r)Y l c (x) and the Hermite/Laguerre basis P\(x). To 
determine precisely which functions Yp{x) appear in the decompositions ( |3. 13a ) 
and ( 3.1 3b| ) , we pick the leading hom ogeneous parts on both sides of the equation. It 
is known from the work of Lassalle |La2 , La3| that the highest-order homogeneous 
part of the multivariable Hermite (type A) and Laguerre (type B) polynomials 
are (with our normalization monic) Jack polynomials with parameter a = 1/go 
in x\,... ,x n and x\ : . . . ,x„, respectively. So, we get upon taking the leading 



POLYNOMIALS RELATED TO THE CALOGERO SYSTEM 



17 



homogeneous part 



|A|/2] 



(3.14a) Jx(x;l/g ) = ]T r 2m Y^_ 2m (x), 



Y 



\\-2m 



e n 



\\-2m 



(3.14b) J A (a; 2 ;l/ 5 o) = 



m— 
|A| 

J2 r 2m Y^_ m (x), 

m—0 



Y 



e n 



B 

|A|— TO' 



where J\(x; l/<?o) and J\(x ; l/<?o) denote the monic Jack polynomial va.x\,. 



and xf, . . . , x\ with parameter a — 1/go [ pt| , M4|. In [Dul], Dunkl provides inver- 
sion formulas for decompositions of the form (B.14a), ( p. 14b ) with which one can 
express the functions Y^_ 2m and m in terms of the homogeneous symmetric 
polynomial in the l.h.s. In our case these inversion formulas become 

(3.15a) 

(3.15b) 

with 



^|A|-2m( X ) 
Y[\\-m( X ) 



TT A 

11 |A|-2m 



J x {x; 1/go), 



TT A 

11 |A|-2m 



rpA 

1 k 



1 



4 m m! [n/2 + d A + |A| - 2m], 
r 2j 



[fc/2] 

E 



^ 4Jj! [-n/2-d A -k + 2]j 



rpA ( t A\m 

~ 1 \\\-2m\ Ij ) 

-(L A y , 



d A = gon(n — l)/2 and 



n 



B 

| A | — m 



J fe — 



1 



4 m m! [n/2 + d B + 2|A| - 2m], 

r 2j 



E iJJ] [_ n /2 - d B - 2k + 2jj 



(L B y, 



d = gon(n — 1) + ng\. Formulas ( 3.13a ), ( 3.13b ) combined with ( 3.15a ), ( 3.15b ) 



render the decomposition of the multivariable Hermite and Laguerre polynomials 
in terms of Dunkl's generalized spherical harmonics in a closed form. 

3.4. Notes, i. Orthogonality. The orthogonality of the multivariable Hermite and 



Laguerre polynomials was stated in [M3 and [La2, La3 . It was pro ved by Macdon- 
ald for the cases go = 1/2, 1 and 2. Recently, Baker and Forrester [ |BF| proposed a 
proof valid for general parameters that exploits the fact that the pol ynomials may 
be seen as limiting cases of the multivariable Jacobi polynomials [0, |Dc2 , BO, |lj. 
In essence this approach should boil down to extending the transition in [Ml] from 
the Selberg integral (whose integrand is th e wei gh t fun ction for the Jacobi polyno- 
mials) to the Mehta-Macdonald integrals ( 2.7a ), ( 2.7b| ) (whose integrands are the 
weight functions for the multivariable Hermite and Laguerre polynomials) so as to 
include also the polynomials of higher degree. 

ii. Norm formulas. Norm formulas for the multivariable Hermite and Laguerre 
polynomials can be found in [M3| and [La2, La3], with again proofs given by Mac- 
donald for go = 1/2, 1 and 2. Baker and Forrester provide a proof of these norm 
formulas for general parameters that hinges on a generating function approach 
(and uses also the orthogonality) [BF|. The expressions for the norms given in 
|M3|, La2, La3 : BF] are written in terms of the Mehta-Macdonald integrals (2.7a), 
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( [2.7b ) and evaluations of Jack symmetric functions at the identity. It can be verified 
with the aid of known evaluation formulas for the Jack symmetric functions due 
to Stanley (see flStJ M4|) that our norm formulas in Theorem 2.2 are in agreement 
with those in [M3, La2, La3, 



BF . 



Hi. Differential equations. The second order differential equation for the mul- 
tivariable Hermite and Laguerre polynomials of the form D^p^ = E^(X)p^ was 
already given by Macdonald [M3] and Lassalle [La2, La3|. In |BF[ a procedure 
is described to obtain the complete system of n-independent differential equations 
starting from the explicitly known differential equations for the Jack symmetric 
functions found by Sekiguchi |@] and Debiard [Del] (see also [M2, |M4[ ). Another 
approach to arrive at such a system of differential equations for p^ is to employ 
Dunkl's differential- reflection operators; see [Ka| for a treatment along this lines 
of the Hermite case. Yet another strategy would be to analyze the limit behavior 
of the system of differential equations for the Jacobi polynomials due to Debiard 
|Dc2 



with respect to the transitions 'Jacobi' — * 'Hermite' and 'Jacobi' 



'La- 
guerre'. All these methods have in common with our approach in Section 2.2 that 
it seems a priori difficult (from a computational point of view) to extract further 
explicit information pertaining to a possible closed form for the higher-order dif- 
ferential equations at the confluent hypergeometric level (cf. also the comments in 
the last paragraph of Section ^2] ). 

iv. Recurrence formulas. Recurrence relations of the type given by Theorem 2.7 
(i.e. the simplest ones, corresponding to the first elementary symmetric function) 
were recently derived independently by Baker and Forrester using a generating func- 
tion for the polynomials [BF]. As it stands, their recurrence formulas are a little 
less explicit than those obtained here since the coefficients are written in terms of 
certain implicitly defined generalized binomial coefficients and furthermore contain 
evaluations of the Jack symmetric functions at the identity. In order make their 
formulas fully explicit (so as to compare with Theorem 2.7) one again needs Stan- 
ley's expressions for the Jack symmetric functions at the identity in combination 
with an explicit representation for the specific binomial coefficients at hand that 
can be found in [ Lai ] . 

v. Rodrigues formulas. Recently, Ujino and Wadati derived Rodrigues type 
formulas for the multivariable Hermite polynomials [UW| following an approach 
due to Lapointe and Vinct who obtained similar Rodrigues formulas for the Jack 
symmetric functions [LV1, LV2|. Such Rodrigues formulas are particularly useful 
when trying to answer questions regarding the structure of the coefficients c\ ^ that 
appear in the expansion of the polynomials in terms of monomial symmetric func- 
tions. For instance, the Rodrigues formulas allowed Lapointe and Vinet to prove 
a weak form of the Macdonald-Stanley conjecture saying that (in an appropriate 
normalization) the expansion coefficients for the Jack symmetric functions in terms 
of monomial symmetric functions are polynomials in the parameters with integer 
coefficients [ LV1 1 . (See [^t], M4] for the Macdonald-Stanley conjecture and vari- 



ous related conjecture s.) A similar statement also holds true for the multivariable 
Hermite polynomials UW]. 
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4. Proofs 

In this section the properties of the multivariable confluent hypergeometric fam- 
ilies stated in Section || are proven by viewing the polynomials as degenerate (lim- 
iting) cases of the multivariable hypergeometric continuous Hahn families (type A) 
and Wilson families (type B) that were investigated in [D3, D4|. 



4.1. Orthogonality properties. In D3, D4| multivariable continuous Hahn and 
Wilson type polynomials were considered that are associated to the weight functions 



(4.1a) A cH (x) = [] 

l<j<k<n 

and 



T(g + i(xj - x k )) 



F(i(xj - x k )) 



[] \T(a + i Xj )r(b- 

l<j<n 



(4.1b) A 



w / 



n 

l<j<k<n 



T(g Q + i(xj - x k )) T(g + i(xj + x k )) 



>< n 

l<j<n 



T(i(xj - Xk)) T(i(xj + Xk)) 
T(a + ixj) T(b + ixj) T(c + ixj) T(d + ixj) 



T(2i Xj ) 



(with go > and Re(a, b, c, d) > 0). Specifically, the multivariable continuous Hahn 
polynomials are defined by the conditions A.l, A.2 in Section 2.1 with A cH (x) (4.1a) 
replacing the weight function A A (x) (1.1a). Similarly, the multivariable Wilson 
polynomials are defined by the conditions B.l, B.2 in Section 2.1 with A B (x) 
( [Tib] ) being replaced by A w (a;) fl4.1bp . 
If we rescale the variables by substituting 



(4.2) 



Xj/P, 



3 = 1, 



and simultaneously perform a reparametrization of the form 



(4.3) 



{P 2 w) 



(p 2 ™7 



gi, 



g'i 



1/2 



(with w, w' > 0, gi,g^ > and P real), then the weight functions A cH (x) (4.1a) 
and A w (ie) ( gig ) pass (upon multiplication by the overall normalization constants 
V A {(3) and V^ifi)) over into 



(4.4a) A A (x) 



V A (P) n 

l<j<k<n 



r(ffo + i/3 l {xj - x k )) 



>< n 

l<j<n 



Tiip-^Xj-Xk)) 
I . 1 



r(^ + i/r^)r( 



m'f3 2 



+ if3 1 x j ) 



and 



A*(x) 



V B (f3) 



n 

l<j<k<n 



r(5o + «/3 - x k )) T(g Q + i/3 1 (x j +x h )) 



(4.4b) 



>< n 

l<j<n 



T(ip-i( Xj - x k )) Tiip-^xj + x k )) 

r(gi + ip- x Xj ) r(gi + 1/2 + ip-ixj) 



T(ip- 1 x j )T(l/2 + iP' 1 x 3 ) 

x r(-^ + ip- x Xj ) r(^— + ip~ x Xj ) 



w'P 
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respectively. The normalization constants V A (f3) and T> B ((3) are introduced so as 
to ensure finite limiting behavior of the weight functions for j3 — > 0. It is not so 
difficult to check (see appendix) — using Stirling's formula for the asymptotics of 
T(z) for \z\ — > oo (see e.g. [ AbSt ]) — that if one takes 

V A {f3) = \P\ 9on ^ n - 1 '>S(zu,l3f n 8{w',P) 2n , 

X> B (/3) = \f3\ 29on{n -^ +2n{si+s '^5(m,f3) 2n 5(w\l3) 2n 



where 



(4.5) 



5{a,p) 



,(l+log(/3 2 a))(l/(/3 2 



-1/2) 



then the weight functions A A (x) ( 4.4a ) and A? (a;) ( 4.4b ) converge pointwise to 



A A (x) (1.1a) and A B (1.1b) for [3 — » provided the parameters of A c (x) are 



related to those of A?(x) by 



(4.6) lu = m + w' and gi = gi 

For our purposes, however, pointwise convergence is not sufficient and we need a 
somewhat stronger convergence result stating that the corresponding measures pass 
over into each other: 



(4.7) 



lim 

13^0 



p(x) A (x) dx\ ■ ■ ■ dx 7l 



p(x) Aa(x) dxi 



(C = A or B) 



where p(x) denotes an arbitrary polynomial in the variables xi, . . . , x n . A proof of 
the limit formula (4.7) can be found in the appendix at the end of the paper. 
Now, let {p A p}\eA an d {P\p}\eA be the bases determined by the conditions 

with the weight functions A A (x) (1.1a) and 



A.l, A.2 and B.l, B.2 in Section 



2.1 



A B ( Lib ) being replaced by A A (x) ( 4.4a| ) and A|(x) ( 1.4b| ), respectively. So, 
up to scaling and reparametrization, the polynomials p A ^(x) amount to the mul- 
tivariable continuous Hahn polynomials (multiplied by /3' A ') and the polynomials 
(x) amount to the multivariable Wilson polynomials (multiplied by /3 2 ' A I) from 



D3 , D4|. It then follows from the defining propert ies for the polynomials (of the 



form A.l, A.2 and B.l, B.2) and the limit formula fl4.7| ) that 

Px(x) - '™ - c 



(4.8) 

and that 

(4.9) (p£,p°)c 



lim 



PxA x ) P°A X ) A^(x)dx 



■ dXr 



(of co urse again assum ing a correspondence between the paramet ers i n accordance 



with (4.6)). Theorem 2.1 is now immediate from limit formula (4.£) and the or- 



thogonality of the multivariable continuous Hahn and Wilson families [D4] (which 



translates into the orthogonality of the bases {p A ^}AeA and {pf^jAeA with respect 



to the weig ht functions A A (x) flOaj ) and A|(x) flOtj)). 

The proof of Theorem 2.2 is based on explicit expressions for the quotients 
(Pa H 'Pa H )ch/(1, 1)cH and (p a V jPa V )w/(1, 1)w — i- e - the ratios of the squared norm 
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of the multivariable continuous Hahn resp. Wilson polynomials and the unit poly- 
nomial (with respect to the L 2 inner product over 1" with weight function A cH (a;) 
resp. A w (x)) — that were computed in \D4\. We have 



/cH 



cH 



= 2_4|A| n 

l<j<k<n 



[go + pf + pf, 1-90+ pf + P c k H }x J+ x k 
[pf + pfA + pf + pf]x j+ x k 



[go + pf-pfA-go 



„cH 



..cHl 

Pk JA,-A fc 



r~cH ~cH i i ~cH ^cHl 

[Pj ~ Pk > 1+ Pj ~ Pk \X 3 -X k 



>< n 

l<j<n 



pf,b cK 



Pf,l 



Pf,l 



hcii 



„cHl 



with pf 1 = (n - j)g + a + b - 1/2, a cIi 
we have 



[pfA + pfh 

a + b - 1/2 and S cH 



6+1/2; and 



(1,1) 



\v 



n 

l<j<k<n 



^0 + Pj + Pfc , 1 - .90 + Pj + Pk]Xj+X k 



[Pj + Pfc> i + Pj + PfcM 3 +A fc 

[5o +pj - P k > 1 - go + Pj - P fe J \j -x k 



pr,i+pr-^k-A fc 



>< n 

i<j<)i 

X 



py ,b w + Pj w ,a w +pf,d w 



[2pf]2A ; 



[i-a w + P f,i-6 w + )0 f ! i-c w + pj v ,i-d 



P j J\f 



[1 + 2pf] 2 A J 



with pf = (n - j)g + (a + b + c + d- l)/2, a w = (a-fe + c- d+ l)/2, 6 W = 



(-a + b + c-d + l)/2, c w = (a + 6 + c + d- l)/2 and 
Here we have adopted the notation 



(-a-b + c + d + l)/2. 



(4.10) 



[ai\i ■ ■ ■ [a p \i 



([a]o = 0, [a]i = a(a + 1) • • • (a + 1 — 1)) for the Pochhammer symbols (and as usual 
|A| = Ai + • • • + A n ). Scaling of the variables and reparametrization in accordance 
with (4.2) and (4.3) leads us to the corresponding expressions for pj^ and Px p, 
which entail for (3 — > (using (O) and (fO|)) 



(4.11a) 



and 



(4.11b) 



(1,1) A 



(^-w n 

l<j<fe<n 



P - j + l)g , 1 + (fc - j - l)ffo]Aj-A t 
[{k ~ j)go, 1 + (fc - j^o^-A* 



x II [! + (»- J)S0]aj 

l<j<n 



(1,1)b 



= (2w) 



-2|A 



n 



[(fc - j + l)g , 1 + (fc - j - Vgohj-x* 

[(fc - j)go, i + (fc - i)5o]A 3 -A k 



l<j<k<n 



l<j<n 
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Combination of the expressions for the ratios in ( i.llaj) an d (|4.11b ) with the 
Mehta-Macdonald formulas for (1, 1)a and (1, 1)b (cf. (2.7a), (2.7b)), produces 
evaluation formulas for (p A ,p A )A and (pf,pf)B that can be cast in the form 
given by Theorem 2.2. Indeed, one easily checks that the norm formulas in Theo- 
rem ^2] are in agreement with the ratio formulas ( 4.11a| ), (4.11b) (using the relation 
T(a+l)/r(a) = [a]i) and, furthermore, that they boil down to the Mehta-Macdonald 
evaluation formulas (2.7a), ( 2.7b| ) for A = 0. In order to verify the reduction to 
( p/Taj ) and ( pJb| ) for A = 0, one uses the identity 

r((* - 3 + 1)30) r(i + (k - j - i) go ) r(i + ngo ) 



ui n 

l<j<fc<n 



2 = 

r((fc - j)go) r(i + (k- j)go) (r(i + go )) n 



(which is derived by canceling common factors in the numerator and denominator 
of the l.h.s. and some further manipulations involving the standard shift property 
for the gamma function z T(z) — T(z + 1)). 



4.2. Differential equations. In [D4| systems of difference equations for the mul- 



tivariable cont inuou s Hahn and Wilson polynomials associated to the weight func- 
tions A cIi (x) ( |0]] ) and A w {x) d4.1b|) were presented. If we rescale the variables 
and reparametrize in accordance with ( |4.2|) and ([Oj) , then the difference equations 
in question pass over into difference equations of the form 



(4.12) 



D%plp = E^{\)p\ 



c 

XP 



(r = l,... ,n; C = A,B), 



for the polynomials P\ a{x) and pf g{x) associated to the weight functions A A (a 



(4.4a) and A? (2;) (4.4b). Here the difference operators D A p 
are the same as in Section 2/2 but with w A (z) and w B (z) replaced by 

(4.13a) 

and 

iz 



and Dfp in the l.h.s. 



(4.13b) 



w B (z) 



w A {z) = (1 + i(3mz){l + ifiw'z) 



1 + 



1 



{iz + (3/2) 

The corresponding eigenvalues in the r.h.s. are given by 



){l + 0wz){l + i0w'z). 



(4.14a) £y^g (A) = (/3 4 ^m'YE r ((p A + Ar) 2 , . 
(4.14b) J B r B / g(A) = (4/3W) r £7 r ((p?+Ai) a , 
with 

(4.15) E r (d,... ,(n,Vr,--- ,Vn) = 

(-ir |J| IIo 



,{p A + X n ) 2 ;{p A ) 2 ,. 
,(^ + A „) 2 ;(pB ) 2 ) 



E 



E 



Vll ■■■f]lr 



JC{1,... ,n} 
0<\ J\<r 



jeJ r<l 1 <---<l r _\ J \<n 



and 
(4.16a) 

(4.16b) 



Pf 



{n- ] )g Q + [3- 2 (- + —\-l/2 1 
(n-j)9o + (go+g'o)/2 + /3- 2 ^ 



1 

2^' 



-1/4. 



In order to derive the differential equations for the multivariable Hermite and 
Laguerre polynomials p A and pf, we shall analyze the behavior of the difference 
equations for p A p and p B ^ when the step size parameter (3 tends to zero. Essential 
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is the behavior of the eigenvalues E^ /3 (X) ( 4.14a ). ( 4.14b ), which is governed by the 
following lemma. 



Lemma 4.1. One has 



(r = l,...,n;C = A,B), 



where E^(X) is given by the expression in Theorem 2.S with to = w + zu' . 



Proof. The proof hinges on the following decomposition of E r (- ■ ■ ; • • • ) (4.15) (cf. 
IM Lemma B.2]) 



E r (Cl, ... , Cn^r, ■ ■ ■ ,Vn) = 

(Cn-Vn)><( E (-l) r " MJ| IIC 



JC{1 n-l} 

0<|,7|<r-l 



j€LJ r < / 1 < • • • <l r — i _ I j I < n 



E (-ir |J| IIo E *x 

jeJ r<ii <•-•<!, | j\<n— 1 



JC{1,... ,«-i} 
0<|J|<r 



with the convention that the sum in the second line equals one if r = 1 and that 
the sum in the third line equals zero if r = n. Notice that these sums are of the 



same type as in the r.h.s. of (4.15) but with one (^-variable less (there is no longer 
dependence on £„). With the aid of this formula and induction on the number of 
variables it is not difficult to infer that 



l l ™P 2r {zu + zu') r M r 



= A„ x 



JC{1,... ,n-l} jGJ 

|J|=r-l 



e nM+ e 



JC{1,... ,n-l}jeJ 

\J\=r 



= e 

jc{i,... .^yjeJ 

\J\=r 

where M = 2 for C = A and M = 4 for C = B. 



□ 



It is immediate from Lemma 4.1 and the limit formula (4.8) that fcWPifl = 



/3 2r E?(X)p% + o(/3 2r ) (assuming of course the usual identification of the param- 
eters via u> — vu + vj 1 and g\ = gi + g[). If we plug this asymptotics into the 
eigenvalue equation ( 4.12Q , then we arrive at the following asymptotic behavior of 
the corresponding l.h.s. for (3 — > 

(4.17) {D%pl p ){x)=^ {D c p c ){x) +o{p 2r )t 

with a certain partial differential operator. In other words, for j3 — > the differ- 
ence equation for p^ „ passes over (cf. Corollary |2.4|) into a differential equation for 



E^(X)p^, where is a certain differential 



the polynomial p'x of the form D^p^ 
operator and E^(X) is of the form giv en in Theorem 2.3 with uj = zu + zu' . 

To complete the proof of Theorem j2.3| it remains to show that the differential 
operator is indeed of the form described by the first part of the theorem (the 
commut at i vity of Df , . . . , £)„ is then an immediate consequence of the commuta- 



tivity (see [|D4 ] ) of D^p, . . . ,D^p). More precisely, we should demonstrate that 
is the leading part of the difference operator D^a in the formal expansion in (3 
and also that its leading symbol is given by the rth elementary symmetric function 
in the partials —d 2 /dx\ 1 . . . , —d 2 jdx 2 n . 
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To this end we first observe that we may replace p^ g in the l.h.s. of (4.17) 



by p^ (using once more that lim^o^A a = f>x)' Moreover, since the polynomials 



{p^}\eA constitute a basis for the space of (even) symmetric polynomials, we have 



in fact that {D^ gp)(x) = (3 2r (D^p)(x) + 0((3 2r+1 ) for arbitrary (even) symmetric 
polynomial p{x) in the variables x\ 1 . . . ,x n . But then the same holds true for 
arbitrary analytic (not necessarily symmetric or even) function of x\, . . . ,x n and 
we have (formally) 

(4.18) D% f) = p» r D? +0{f3 2r+1 ). 

Here we have used the fundamental property that the vanishing of a (linear) partial 
differential operator on the space of (even) symmetric polynomials implies it be zero 



on an arbitrary (analytic) function (i.e., all its coefficients must be zero). In | D1 
Appendix C] a proof of this general property was given in a trigonometric context. 
(Specifically, we assumed there that the differential operator vanishes on the space 
of symmetric polynomials in sin 2 (cci), . . . , sin 2 (x n ).) The present case follows after 
an appropriate substitution of the variables turning the relevant space of trigono- 
metric polynomials into the space of (even) symmetric polynomials in x\, . . . ,x n . 
(Specifically, the type A case is recovered by the substitution sin (xj) — ► Xj and 
the type B case by sin 2 (xj) — > x 2 .) A consequence of this general property is that 
the leading part of the formal expansion of the difference operator D^g in (3 is 
completely determined by its action on the (even) symmetric polynomials. (This 
excludes the (a priori) possibility of a lower-order leading part in the formal expan- 



sion (4.18) corresponding to a term determined by a nontrivial differential operator 
that vanishes on the space of (even) symmetric polynomials.) Notice also that this 
property may be used to arrive at an alternative proof for the commutativity of 
Dp, . . . , D„. Indeed, the commutators [D^, D^] obviously vanish on the simulta- 
neous eigenbasis {p a }agA (and hence on all (even) symmetric polynomials), from 
which it then follows that they be zero identically. 

To determine the highest-order symbol of the leading differential operator in 
( 4.18] ) , we use that the functions v c , w c governing the coefficients of the difference 



operator D^ g are of the form 1 + 0({3) and that for v c ,w c = 1 



r, (v c w c = l) x -i-r / £-5- -S.-0- 

(4.19) D®p (v ,w = 1} ]T II ' 

JC{1,.-. ,n}jeJ ^ J 

\J\=r 

= /3 2r ((-D r E Il^r) +^ 2 " +1 )- 



jc{i,... ,»>ieJ 3 

\J\=r 



(Notice in this connection that in the situation where v c and w c are taken to be 



equal to one, the operators D^g (2.8a), (2.8b) reduce to 



,A _ V- ( _ 2) r-\J+\-\J-\ ( n ~\J+\ ~ \ J -\ \ A 

r - l3 ~ J+ ,_h ( } [r-\J + \-\J-\) e 

J + HJ_=0, |J + | + |J_|<r 
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and 



U r,@ — 



(-2) 



r-\J\ 



JC{1 «}.0<|.7|<r 



\J\ 

\J\ 



which both can be rewritten in the form given by the first line of (4.1S ).) It then 



follows (i.e 
part Z)J? in ( [4.1 



from the asymptotics v c , 
is of the form 



0{/3) and ( ff,19| )) that the leading 



D? = (-iy E 



JC{1,... , 
|J|=r 



as advertised. 

Thus far we have shown that the statements of Theorem |2.3| h old when taking 



difference operators D^ g (2.8a), ( [2.8b ) with w c (z) given by (4.13a), ( 4.13b| ), where 



corresponds to 



w + w — uj and gi + g[ = g\ . The formulation in Theorem |2.. 
choosing the specialization w = u>, vj' = and gi = gi, g[ = 0. 

The symmetry of the differential operator D^ Q with respect to the inner product 
(•, -)c in the space of permutation-invariant (type A) or permutation-invariant and 
even (type B) polynomials (Corollary 2.5) stems from the fact that the operator in 
question is diagonal on an orthogonal basis (viz. {p^}asa) with eigenvalues that 
are real. 



4.3. Pieri type recurrence formulas. In [D4] Pieri formulas for the multivari- 
able con tinuo us Hahn and Wils on polynomials associated to the weight functions 
A cH (x) ( 4.1a| ) and A w ( x) (|4.1b ) were introduced. After rescaling and reparametriz- 
ing in accordance with ( |4.2[ ), ( |4.3| ) , one arrives at the corresponding Pieri formulas 
for the polynomials p A a and p\ s (which — recall — are determined by the conditions 



A.l, A. 2 and B.l, B.2 of Section 2.1 with the weight func tions A A (x) ( 1.1a ) and 
A A (x) ( |l.lb| ) replaced by A A (x) ( |4.4aD and A|(:c) ( |4.4bD ). In the simplest case 
the resulting Pieri formula takes the form 



(4.20) EC 8 (x)P° 3 (x) 



l<j<n 

E V c J .p(p c + \)(pC_ e] ^x)~pC a (x)) 



l<j<7l 

\-ej SA 



with 



\<k<n.k=£j 



The functions v c , w c are given by 



v A (z) = 1 + 2°, ^ A (z) 



{a A + z){b A + z) 



I 



f ,B( z ) 1 + 50 B( , (a B + zW + z)(c B + z)(cP + z) 

[ > ~ + z ' [ >~ 2z(l + 2z) 
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with 



/r 2 (- + — )-iA 

ZD ZD 



/r 2 (- 



1/2 



and 



-B 



1 



1 



= ( gl + g' 1 )/2 + /r 2 (- + -)/2-i/4, 

6 B = ( gl + gi)/2 - /T 2 (- + i 7 )/2 + 3/4, 
zd w 

c B = (gi - gi)/2 + /3' 2 (- ~ A)/ 2 + 1/4, 

1 1 



gi)/2 -/?-'( 7 )/2 + l/4. 

w zd 



The vector p c has components given by (4.16a), ( 4.16bj ) and the multiplyer in the 
l.h.s. of (4.20) is given by 



(4.21a) 

(4.21b) 

with 
(4.22) 



l<j<n 



l<j<n 



9 



■pf) 



(pf) 2 ) 



pf 



pf 



(n - j)go + gi 



(n - j)go + l/{P 2 w) and 
In the Pieri formula we have furthermore employed the normalization 
(4.23) 
with 



c A,/3 



n 



l<j<n L" 1 
~A j_ ~A 



1 , k<n t- 90 + P f + Pkhi+^k [90 + pf - Pk]\j-\ h 



= (-i//? 2 ) |ai n 



2A, 



l<j<n 



n 



+ P J B ,6 + p j B ,C + p j B , d + Pj B ] Aj 



l<j<fe<« t- 90 + Pf + PkU+^ [-90 + Pf - PfcU 3 -A fc ' 

The recurrence relations of Theorem 2.7 can be recovered from ( |4.2C ) for j3 
To see this, one first observes that 



(4.24) \jm(-£-) W P£ fi {x)=P£(x) ) 



YimP* (x)=P%(x) 



with Px(x) given by (|2.12|) . These limits follow from (fl^) and the fact that 



(4.25) 



' JM A - A lira l! ''' 



z \,l3 — c \ 



with c A and c B given by (|2.13aD and ( |2.13b| ). (As usual, we assume an identification 
of the parameters of the form lu — zd + w' and g\ = gi + gi-) For type A, 
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multiplication of ( 4.2C ) by i/3(^)' A ' leads for (3 — > to the first recurrence relation 



of Theorem 2.7 The second (i.e. type B) recurrence relation of Th eorem 2.7 is 
obtained similarly by sending /3 to zero in the type B version of (4.20) after having 
multiplied both sides by the factor [w + w')(3 2 . 

To derive these limit transitions for the recurrence relations we have used that 
for (3 -> 

i°(pf +Pk+ A, + Xk),i> c (~pf - p£ - Aj - A fc ) = 1 + 0(/3 2 ), 



go 



( fc - j)9o + Aj - A fe 



and 



™ (Pi +Aj) 



/3 2 7Z7 



(l + 0(/3 2 )), 

- j)go + Aj 
2(tj7 + m') 



(l + 0(/3 2 )), 



- B (P, B + A,) = ("-^+8i+8i + 1 /2 + A 3 (1 + Q 2 



.u x ,-^-^ + A i ( i + (/3 2 )). 



! (-pf-A,) = 



/3 2 (tu + ro') 

Notice to this end also that in the case of type A, a divergent term in the l.h.s. 
of the recurrence relation (4.20) originating from the factor — . p A (cf. (4.21a)) 
cancels against a corresponding divergent term in the r.h.s. originating from the 
factor in front of P\p of the form — • V J A 3 (p A + A). (That the divergent terms on 
both sides indeed cancel is seen using the identity J^j Tikjtji^ + 9o/(Cj — 00) = n -) 
In general the recurrence relations for the polynomials P\ p ( |4.23| ) induced by 
|D4|1 become 



(4.26) E^(x)P^(x) 



E u%, r _ vn { P c + \)V^Ap c + \)pv +eEj ^{x), 



JC{1,... ,»}, 0< J|<r- 

£j=±l, A+e eJ £A 



, n, with 



E 



fee {+1,-1}), 



i6J 



j.j'ej 
j<j' 



n« c fe-Ci+Cfe)« c (£iO-a), 



^(0 = 

(-1)" 



E (n* C ( £ ^0 II « C (ei0+ei'Ci')*°(-e«0-eKj'-l) 
I] * C (eiCi + Cfc)« C (ejC«-C*)) 



icii, |i|=p i£L 

ej=±i, zeL 



k£K\L 
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and 



(4.27a) E A {x 



) = (-1) 



(4.27b) E*Jx) = 



E nf 

JC{1 n} M r<Z 

0<|J|<r 



<— <2r- 



••Pl r 



E 

jc{i,... ,»} je./ 

0<\J\<r 



TI- 
LL 02 



E 

r<£i — I ji <t 



For r = 1 the recurrence formula in ( 4.26| ) specializes to that of ( 4.20| ). 

It is not difficult to see that the recurrence relations for the multivariable La- 
guerre polynomials characterized by Theorem 2.S and Theorem |2.9| f ollow from the 
type B version of (4.26) for j3 — > 0. Indeed, multiplication of ( 4.26| ) by the factor 
P 2r {va + vj'Y and sending j3 to zero readily leads to the Laguerre type recurrence 
relations. The verification of this assertion hinges on the second limit formula of 

3 displayed above, and the fact that 



(4.24), the asymptotics for v 



lim B 2r E, 
p-> o 



(-If 



E 

JC{1,... ,n}jeJ 
\J\=r 



n 



For type A the transition 8 — > is substantially more complicated due to the 



singular nature of the terms in ( 4.26 ). Specifically, the multiplyer E A 3 {x) ( 4.27a ) 



consists of a linear combination of the elementary symmetric functions in x\, . . . ,x n 
up to degree r. The coefficients in this linear combination have a pole at (3 = 
0, the order of which is reversely proportional to the degree of the elementary 
symmetric function in question (notice that pf — 0(/3~ 2 )). Hence, for f) — * 
the contributions of the lower-degree elementary symmetric functions to E^(x) 
become predominant. To get rid of these lower-degree divergent terms, we take 
an appropriate linear combination of the recurrence relations (4.26) that cast them 
into a system of the form 



(4.28) 



with 



E 

./e{i,... , 
\J\=r 



II*: 



E 



W A P A 

vv eJ, J c ;r,p r X+e e 



(x) 



JC{1,... ,n},0<|J|<r 

,=±1, e e j+A£A 



Basically, this boils down to passing from Picri formulas corresponding to the sym- 



metric functions E A Jx) ( 4.27a ) to Pieri formulas corresponding to the elementary 



symmetric functions Y2\ j\= r TijeJ x i ^ subtracting from the rth Pieri formula in 



( |4.26| ) a suitable linear combination of the Pieri formulas corresponding to E A g (x) 
with s < r (and multiplication by an overall factor). The coefficients of the terms 
in the r.h.s. of ( |4.26| ) labeled by index sets J with | J\ = r are invariant with respect 
to such changes in the l.h.s. (up to an overall factor (ij3) r and a factor caused by 
the change of the normalization P A g — > P\g)- More precisely, we obtain that for 

\J\ 



r the coefficient W£j j c . r g in the r.h.s. 



of (4.2 



W. 



eJ, J c ;r,l3 



eJ, J 



The type A version of Theorem |2.8| then follows for /3 
(4.24) and the above asymptotics for v A and w A ). 



is given by 

(P A + A). 
(using the limit formula 
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It remains to verify the normalization properties of Pjp(x) stated in Theorem 2.6 



These properties are a consequence of the fact that P^(i(3p c ) = 1 (see the remarks 



in [D4, Sec 6]). Specifically, by sending f3 to zero in the relation Pjp JififP) = 1 



we arrive at the normalization properties of Theorem 2.6. For type B this is im- 
mediate from the limit formula (4.24); for type A this is seen by noticing that 
lim/3^0 P\ /j(*/3p A ) picks up the highest-degree homogeneous part of P^g(x) evalu- 
ated in x = 1 (here we use that p A = -^kp + 0(1) together with the limit formula 
( 1.25| )), which is equal to lim Q ^oo a _ ' A 'P^(al). 



Appendix: Convergence of the weight functions 



In this appendix it will be shown that the weight functions A^(ir) ( 4.4a ) and 



Ag (x) (4.4b) converge for (3 



f\u>j y ft-^uj / tuuvcigu iui )j —r to the weight functions A (x) (1.1a) and A (x) 
(1.1b), respectively. More precisely, we will prove the somewhat stronger result 
that 



(A.l) 



lim 

p^o 



p(x) A^(x) dxi 



p(x) A c (x) dxi ■ ■ ■ dx n (C = A or B), 



where p{x) denotes an arbitrary polynomial in the variables x±, . . . ,x n . In (A.l) it 
is understood that the parameters of A c (x) ( |l.la| ) , ( 1.1b ) and A^ (x) ( |4.4a| ) , ( |4.4b| ) 
are related via the identification lo = w + w' and g\ = gi + g[ . 

Let us start by inferring the pointwise convergence of the weight functions. For 
this purpose we use the limit formula 

1 



(A.2 
with 



lim l <$(<*,/?) |r( 
P^o ap z 



+ i/T 1 y)| =exp(-ay 2 /2) 



{a > 0) 



Sfap) = /^ e (i+i°g(^ 2 ))(«- 



and the limit formula 



(A.3) 



lim 

p^o 



a Y{a + b + i(3- l y) 



r(6 + z/3- 1 y) 



M>o), 



where in both formulas it is assumed that y and (3 are real. By applying (|A. 2|) 



and (A.3) to the factors of A-g(x) (4.4a) and A?(x) (4.4b), one readily sees that 



f or (3 — > these weight functions converge pointwise to A (x) ( 1.1a ) and A 



( 1.1b) as indicated. The normalization factors of the form 5(a,f3) and \(3\ a in (A.2) 
and (A.3) ensure a finite and nontrivial limit; the factors in question have been 



collected in the weight functions A^(x) ( 4.4a ) and A^(x) ( 4.4b ) into the over 



all normalization constants V A (/3) = \/3\ n ( n - 1 '> 90 6(zu, (3) 2n 6(zu', f3) 2n and V B {fi) = 

| /3 |2n(n-l) S0 +2n(gi+gi) ( J( W) fffn 

The limit formulas ( A.2 ) and (|A.3| ) may be verified with the aid of Stirling's 
formula for the asymptotics of the gamma function for large values of the argument, 
which reads (see e.g. AbSt , |oi|l ) 

(A.4) 



T(z) = (27T) 1 / 2 e^z 2 " 1 / 2 • exp(i?(z)) 



:!() 
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with R(z) = 
z = a^p- 2 

(A.5) 



0(l/|z|) f or \z \ — > oo in the sector | arg(z) | < n. Substitution of 
- i(3~ 1 y in ( A.4 ) entails that for {3 — > 

1 



5(a,0) 



(l + o^V) 21 " 3 



e ' 3 



which implies (A.2). In a similar way one concludes from (|A 



rctan(Q/3y) (l + 0(/3 2 )), 

that for (3 -> 



(A.6) 



r(& + i/?- 1 ^) 



which implies (A3). 

Let us next demonstrate that the pointwise convergence of the integrands car- 
ries over to the convergence of the integrals by invoking Lebesgue's dominated 
convergence theorem. For this pur pose it is needed to dominate (the absolute value 
of) the integrand in the l.h.s. of (Al) uniformly in (3 by an integrable function. 
We will do so by bounding individually the factors comprising the weight func- 
tion. Specifically, it turns out that factors of the form 5(w,(3) |T(^2 + i(3~~ 1 Xj)\ 
(and 5(w',f3) |r( L -I- i/3 Xj)\ ) may be dominated by an exponentially decay- 
ing function and that the remaining factors — which consist of ratios of the form 
|/3 9 °F( 5o + if3-\x 3 ± x^/Tiip-^Xj ± a*))], \(3 Sl T( gl + i^x^/T^x^ and 
|/3^T(gi + l/2 + i/3~ 1 Xj)/T(l/2 + if3~ 1 Xj)\ — grow at most polynomially in the vari- 
ables X\, . . . ,x n . Hence, the integrand on the l.h.s. may be dominated by an 
exponentially decaying function and (A.l) follows from the pointwise convergence 
of the integrands by the dominated convergence theorem as advertised. 

To validate the above claims regarding the bounds on the growth of the factors 
constituting the weight functions Ap(x) (4.4a) and Aj^(x) (4.4b), we need a precise 
estimate for the error term R(z) appearing in the Stirling formula joj] 

1 



(A.7) 



\R(z)\ < 



= arg(z)). 



12|z|cos 2 (6>/2) 

(This estimate is valid in the whole sector |arg(z)| < tt.) By substituting z = 
a -1 /3~ 2 + if3~ l y in the Stirling formula (A.4) (where a is assumed to be positive) 
and combining with the error estimate (A.7), we find (cf. QA.5|)) 



(A- 



with 



(A.9) 



S(a,(3) \T{^-+ % ^y)\<e~ F ^Gp, 

a(3 z 



Fp{y) 



y 

— arctan(a/5y) 
P 



i 



2a(3- 



log(l + a 2 /3V) 



and Gp = e Q/3 / 6 . (In our situation cos 2 {9/2) > 1/2 in view of the fact that the 
real part of z = + i(3~ 1 y is positive.) Differentiation of Fp(y) with respect to 
y yields 

(A.10) d y Fp(y) = /T 1 arctan(/3ay), 

which shows that Fp(y) is nonnegative as an increasing/decreasing function for y 
positive/negative with Fp(Q) — 0. From the asymptotics for \y\ — > oo one further- 
more sees that the factor exp(— Fp{y)) decays exponentially. A little more precise 
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analysis reveals that for < [3 < 1 

To obtain the exponential bound on the tail we have used: (i.) that for < (3 < 1 



Fpfl/a) = — arctan(/3) 



1 



log(l + /3 2 ) 



2a(3 2 

> (7r/4-log(V2))/ a > l/(3a), 

(ii.) that for y > 1/a and < (3 < 1 

(») > (fy-fy) (1/a) = /T 1 arctan(/3) > tt/4 > 1/3, 

and (iii.) that Fp{y) is even in y. We conclude from JA.8| ) and ( |A .11 ) that for 
< (3 < 1 the factors in the weight function (x) (C — A, B) of the form 
5(m, f3) \T(-^p2 + if3~ 1 xj)\ (and S(ru', (3) |r(^7^2 + 0~ lx j)\ ) may indeed be domi- 
nated uniformly in (3 by an exponentially decaying function of Xj . 

It remains to check that the rest of the factors — which consist of ratios of gamma 
functions — can be dominated by a function that grows at most polynomially in the 
variables. To this end we shoul d fin d bounds on the gamma function ratios of the 
type appearing in the l.h.s. of (A. 3). Notice that for a 6 N we have 



(A.12) 



a T(a + b + if3- 1 y) 



T(b + if3-iy) 



J[ \(m + b)(3 + iy\, 



m=0 



which is easily dominated uniformly in < [3 < 1 by a function with polynomial 
growth in y (take e.g. the function ((a + b) 2 + y 2 ) a ^ 2 )- The case of general positive 
(not necessarily integer valued) a is a little less straightforward; it will be addressed 
here with the aid of the integral representation 



(A.13) 



T(a + z) 



z exp 



1 



1 



Re(z) > 0, 



which can be obtained by integrating Gauss' integral formula for the psi function 
ij){z) = T'(z)/T(z) [|AbSt|, 01 



t 



1 - e" 

-t 



(A.14) 4>(z) = 

Jo 

(using also that log(z) = J °°t~ 1 (e 
moment assume that b is positive. Then substitution of z = b + if3~ 1 y in (A.12) 
entails 



dt Re(z) > 

tz )dt for Re(z) > 0). Let us for the 



(A.15) 



a r( a + 6 + i /3- 1 y) 



T(b + i(3- x y) 
{(3 2 b 2 



2/ 2 ) a/2 exp 



e~ bt cos(yt//3) a 



1 



The integral within the exponent is bounded by a constant with value 

dt 

T 





1 - e~ at 


/ e- bt 


a t 


/o 


1 - e"* 
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and the factor in front is smaller than (6 2 + y 2 ) a / 2 for < (3 < 1. The case that b 
becomes zero can be reduced to the previous situation with positive b by means of 
the identity 

Tja + ip-iy) = Tja+l + i^y) ( iy 
( ' ' T{i/3-iy) r(l + i/3-iy) \a0 + iy, 

The upshot is that for < /? < 1 the factors \(3 9o T(g a + i/3~ 1 (x J ■ ±x k ))/r(ip- 1 (x j ± 

x k ))\, IpsiT^+ip-^/Tiip^Xj)] and |/3g^^(g' 1 + l/2+^/3- 1 x,)/^(l/2+^/?" 1 x J )l 
can be uniformly dominated in [3 by a function that grows at most polynomially in 
the variables x\ , . . . , x n , which completes the proof of ( A. 1 ) . 
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